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Abstract 

A quantization scheme for the phenomenological Maxwell theory of the full elec- 
tromagnetic field in an inhomogeneous three-dimensional, dispersive and ab- 
sorbing dielectric medium is developed. The classical Maxwell equations with 
spatially varying and Kramers-Kronig consistent permittivity are regarded as 
operator- valued field equations, introducing additional current- and charge- density 
operator fields in order to take into account the noise associated with the dis- 
sipation in the medium. It is shown that the equal-time commutation relations 
between the fundamental electromagnetic fields E and B and the potentials A 
and (p in the Coulomb gauge can be expressed in terms of the Green tensor of the 
classical problem. From the Green tensors for bulk material and an inhomoge- 
neous medium consisting of two bulk dielectrics with a common planar interface 
it is explicitly proven that the well-known equal-time commutation relations of 
QED are preserved. 



1 Introduction 



It is well known that the quantum statistical properties of electromagnetic fields 
including their interactions with atomic systems can be strongly influenced by 
the presence of dielectric bodies. Typical examples are the Casimir effect |]. 0- 
the modification of the spontaneous emission rate of excited atoms || |, H in the 
presence of dielectric media, and the degradation or improvement of nonclassical 
properties of light propagating through optical devices, such as cavities, beam 
splitters, wave guides etc., which typically can be regarded as dielectric bodies|| 
[7|, ||. Therefore it has been of considerable interest to formulate QED on a 
dielectric-matter background. Various quantization schemes have been proposed 
for dispersionless [g |H], |TJ, 0, [Ej, [L|, p|, @ , dispersive [0, [18], [L9|, 0, |lj, g| [EJ , 
and nonlinear |18|, [TJ], 2IJ ^5, |24J, |25], |26| dielectrics. However, most of these 
quantization schemes run into difficulties when an absorbing medium is attempted 
to be included in the concept, which is crucial for studying propagation effects 
and keeping the theory consistent with the causality principle. 

The problem has been considered by a number of authors [^7, |28|, 29, 30, [31 
32], 0, 0, ||, |§. In p 



a fully canonical quantization scheme for the 
macroscopic electromagnetic field in a linear harmonic-oscillator bulk material 
is developed which is based on the Hopfield model of a dielectric P5| . The 
electromagnetic field is coupled to a harmonic-oscillator polarization field that 
interacts with a continuum of harmonic-oscillator reservoir fields. The resulting 
Hamiltonian, which is a bilinear form of bosonic fields, is diagonalized in two 
steps - first the polarization-reservoir part and after that the total Hamiltonian. 
The scheme is much more involved when it allows the electromagnetic field to 
be in an inhomogeneous medium, as it is the case in practice, and much effort 
must be made to perform the diagonalization even for simple dielectric-body 
configurations pTj . 

Another approach to the problem of including losses in the quantization 
scheme is the method of Green function expansion [BO, |4l|, which can be re- 



garded as a natural extension of the familiar method of mode expansion (which 
only applies to strictly nonabsorbing media) to arbitrary Kramers-Kronig consis- 
tent media. The approach, which resembles, in a sense, the method of (operator) 
Langevin forces 0, [31], [52], ^3|, directly starts with the Maxwell equations for 
the macroscopic electromagnetic field, including the dielectric displacement vec- 
tor and a (phenomenologically) given permittivity. The quantization of the ra- 
diation field is based on the classical Green function representation of the vector 
potential, identifying the external sources therein with the noise sources that are 
necessarily associated with the losses in the medium and replacing the c-number 
sources with operator-valued ones such that the equal-time basic commutation 
relations of QED are preserved. The advantage of the method is that the calcula- 
tion of the Green function is - similar to the determination of the mode structure 
in the standard scheme - a purely classical problem. The Green function is essen- 
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tially determined by the permittivity of the medium, which is a space-dependent, 
complex function of frequency. The configuration of the dielectric bodies is de- 
scribed by the dependence on space of the permittivity, and the effects of disper- 
sion and absorption, respectively, are described by its real and imaginary parts. 
It is worth noting that they are not independent of each other, but they must 
satisfy Kramers-Kronig relations, because of causality (see, e.g., |44], [45]]). The 
Green function method has been proved correct for radiation in 3D bulk ma- 
terial [JD, E6[ and in ID multilayer structures [4"0, O, 46], and applications to 



various problems have been studied (e.g., ground-state field fluctuations [41], flTfl, 
photonic wave packets at dielectric barriers nonclassical-light propagation 



in dispersive and absorbing dielectrics [49], p0|). 



The aim of this paper is to extend the Green function method to the quantiza- 
tion of the electromagnetic field in a dispersive and absorbing 3D inhomogeneous 
medium. For this purpose, both the transverse and the longitudinal parts of the 
electromagnetic field vectors must be included in the analysis in an unified man- 
ner. Relating the overall noise current to a bosonic basic field, the electromagnetic 
field operators can be expressed - through the dyadic Green function of the classi- 
cal problem - in terms of this field, and all the fundamental electromagnetic-field 
commutation relations can be expressed in terms of the Green function. Using 
earlier results of the calculation of the classical 3D Green function for multilayer 
systems [pT], |52|, [53], EM , we explicitly consider an inhomogeneous medium that is 



composed of two bulk dielectrics with a common planar interface. 

The paper is organized as follows. In Sec. |2] the quantization scheme is devel- 
oped. The Green function expansions of both the electromagnetic field vectors 
and the vector and scalar potentials in Coulomb gauge are given, and the funda- 
mental commutation relations of QED are studied. In Sec. |3| the theory is applied 
to a bulk dielectric, and it is shown that the equal-time basic commutation re- 
lations of QED are preserved. In particular, earlier results for the transverse 
electromagnetic field are recovered. Quantization of the electromagnetic field in 
an inhomogeneous medium that consists of two bulk dielectrics with a common 
interface is studied in detail in Sec. [4] Finally, a summary and conclusions are 
given in Sec. [5| Lengthy mathematical derivations are outlined in appendices. 

2 Quantization scheme 
2.1 Classical Maxwell equations 

Let us start by writing the phenomenological Maxwell equations in the frequency 
domain as 

V-B(r,w) = 0, (1) 
V-D(r,cu) = 0, (2) 
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VxE(r,w) = icuB(r,uj), (3) 
VxH(r,w) = -iuB(r,u), (4) 

where we have assumed that no "visible" charges and currents are embedded in 
the background medium. The Maxwell equations must be supplemented with 
constitutive relations, which for linear dielectric media are usually given by 

D(r,w) = e e(r,w)E(r,w), (5) 
B{r,u) = /x H(r,u;). (6) 

Here, the (relative) permittivity, which is a complex function of frequency, 

e(r,u;) = e R (r,uj) +iei(r,u), (7) 

is also allowed to be varying with space in order to model inhomogeneous media. 
For causality reasons, the real and imaginary parts of the permittivity, which 
are responsible for dispersion and absorption, respectively, are uniquely related 
to each other through Kramers-Kronig relations, i.e., dispersion and absorption 
are intimately linked. It can be shown that e(r, Q) as function of the complex 
frequency Q is analytic and has no zeros in the upper complex half-plane, and 



e(r, fi) — > 1 if \Q\ — > oo |44. 45|. The fields in the time domain are obtained by 



Fourier transforming the fields in the frequency domain, e.g., 

POO 

E(r,t)= duje' wt E(r,uj)+c.c., (8) 



o 

and B(r,t), D(r,t), and H(r,t) accordingly. 

The Maxwell equations (|I|) - (J|) together with the constitutive relations (|5|) 
and (0) cannot be transferred to quantum theory by simply regarding the elec- 
tromagnetic field vectors as operator-valued quantities, otherwise the operators 
would be damped to zero. This is not surprising, because equations of the form 
given here violate, in general, the dissipation-fluctuation theorem, which states 
that damping is always connected with additional noise. In other words, even 
the classical equations are equations for the field averages but not equations for 
the "naked" fields, and therefore they cannot be used to study the statistics of 
fluctuating fields, such as thermal fields. Hence, transferring the above given 
equations to quantum theory can only yield equations for the (now quantum- 
mechanical) expectation values of the fields (which of course can be damped to 
zero). The noise that is unavoidably associated with absorption can be described 



by introducing a corresponding source term in the Maxwell equations Jl2, 55] , 
which can be thought of as to arise from a noise polarization in the constitutive 
relation between the dielectric displacement vector and the vector of the electric 
field strength (see, e.g., @]). 

D(r,w) = eo e(r,^)E(r,^) + P(r,cu). (9) 

Before specifying the noise source let us first turn to quantum theory. 
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2.2 Quantum Maxwell equations 



The Maxwell equations ([!]) - (§) together with the constitutive relations (f|) and 
(HD can be transferred to quantum theory, regarding the electromagnetic field 
vectors and the noise polarization field vector as operators: 



0. 



(10) 

e e(r, a;) E(r, a;)] = p(r,u), (11) 

VxE(r,w) = iwB(r,w), (12) 

V x B(r, u) = -iujp e e(v, u) E(r, u) + /i j(r, w). (13) 

Here, the operator noise charge density p and the operator noise current density 
j are introduced, which are related to the operator noise polarization P as 

p{r,u) = -V-P(r,^), (14) 
j(r,w) = -iwP(r,w). (15) 

It follows from Eqs. ([14]) and (p~5|) that p and j fulfill the equation of continuity, 

V • j(t,u) = iujp(r,uj). (16) 

The electric-field strength operator E(r) (in the Schrodinger picture) is defined 
in terms of the Fourier transform E(r, u) as 



E(r) 



dujE(r,uj) + H.c, 



(17) 



and similar relations hold for B(r), D(r), and H(r). 

As already mentioned, the source terms p and j are closely related to the noise 
associated with the losses in the medium, which themselves are described by the 
imaginary part of the permittivity. Following jlO], |4"I|, we relate j to a bosonic 
vector field f as 



J(r,wj 



Hoc' V Treo 



h 

— ei(r,w)f(r,u;), 



5ij5(r — r')8(u — lu'), 







(18) 

(19) 
(20) 



Obviously, in the Heisenberg picture the basic operator field evolves as f(r, u,t) 
= f (r, lu, t') exp[— iu(t — t')], which is governed by the Hamiltonian 



H = J d 3 v ^ duhutf(r,u) ■ f(r,w). 



(21) 



The system of equations (|T0| ) - ([13]) together with Eqs. ( |TT| ) - fl2~T|) is complete, 
i.e., further equations are not required. In particular, all the electromagnetic- 
field commutation relations are uniquely determined from the equations given. 



It should be pointed out that - in contrast to J4(J - the current density j is 
not transverse, because the whole electromagnetic field is considered. Hence, the 
vector field f introduced here is not transverse as well, and the spatial 5 function 



in Eq. ( |i9|) is an ordinary 5 function instead of a transverse one. 
2.3 Integral representation of E(r, uS) and B(r, uj) 



Equations fll2|) and (0) imply that the electric field E(r, to) obeys the partial 
differential equation 

uj 2 

V x V x E(r, uj) -e(r, u;)E(r, uj) = iujUojU, uj), (22) 

cr 

whose solution can be represented as 

uj) = iufio J d 3 sG i:j (r,s,uj)l.(s,uj), (23) 

where j is given by Eq. (|T8|), and Gij(r, s, uj) is the dyadic Green function (Green 
tensor) of the classical problem. It satisfies the equation 

,2 



dld r m -5 im lA r + ^e(r,co) 



G m j(r, s, u) 



= <J fi< y(r-s), (24) 

together with appropriate boundary conditions. In particular, it must vanish at 
infinity. The notation <9[ means d/dxi, and A r = <9[<9[ (here and in the following 
the summation convention is used). When the electric field E(r,u) is known, 
then the magnetic field B(r, uj) can be obtained as, on recalling Eq. (|T^), 

B(r, uj) = (iuj)- 1 ^ x E(r, uj). (25) 

In this way, the electromagnetic field operators can be expressed in terms of the 
bosonic field f(r, uj). 

2.4 Commutation relations 

Let us consider the (equal-time) commutation relations between the fundamental 
fields E(r) and B(r). Recalling the definitions of E(r) and B(r) [see Eq. (|17|)1, 
using Eqs. (^) and (|25|) together with Eq. (|18D, and applying the commutation 
relations (|19| ) and (|20D , we find that 

[E i (r),B k (r')] = —e kmj ff m dw 
7reo Jo 

x ^ J d 3 s ei(s, uj)G a (r, s, uj)G* jl (r', s, uj) - c.c, (26) 
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where ekmj is the Levi-Civita tensor. In order to simplify Eq. fl26|), we first note 
that the relation 



OJ 



J d 3 s ei(s, u)G u (r, s, cj)G^(r', s, u) 



1 

27 



^•(ry^) - G*-(r, r',w) 



(27) 



is valid (see App. |A|). Further, from Eq. ([2~3|) and the relation e*(r,u;) =e(r, — co>) 
it follows that 



Combining Eqs. ( |26| ) - (p8|), we derive 

h 

7re 

Similarly, we find that 

~EAr),E k (r') =0 



(28) 



da;— Gij(r,r ,u). 

-oo C 



(29) 



(30) 



which is in full agreement with QED (see, e.g., |57| ). Equation (|29| ) reveals 
that the commutator between the electric and magnetic fields can be expressed 
in terms of a single frequency integral of the Green function multiplied by the 
frequency. In order to calculate this integral, knowledge of the Green function is 
required. Note that a single pole at u = has to be treated as principal value. 

2.5 Vector potential and scalar potential 

2.5.1 Potential equations and integral representations 

It is often necessary to use electromagnetic potentials. In the frequency domain, 
the vector and scalar potentials A and <p, respectively, are related to the fields as 



B(r,a>) = V x A(r,u), 

E(r,c<j) = iuA(r, u>) — Vy?(r, u). 



(31) 
(32) 



Substituting in Eqs. ( |ITD and flT3] ) for the fields the potentials according to 
Eqs. (P) and (||) then yields 



V- [e(r,a;)V0(r,u;)] 



-— + iujV ■ [e(r, oj)A(r, u)], 



(33) 
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V x V x A(r,u) -e(r,u;)A(r,u;) 

er 

itu 

= y"oj(r, w) + -ve r, u)V<p(r, u). (34) 

In Coulomb gauge, 

V ■ A(r,cj) = 0, (35) 

Eq. ( |32"D corresponds - in the sense of the Helmholtz theorem - to a unique 
decomposition of the electric field E into a transverse part iuj A and a longitudinal 
part — V<£>. Hence we may write 

Mv,cu) = (iu)- 1 J rf 3 s^(r-s)4(s,u;), (36) 
d^(r,cu) = -J rf 3 s4(r-s)4(s,^, (37) 

where 5 _L (r) and 5"(r), respectively, are the transverse and longitudinal 5 func- 
tions (see, e.g., [§), 

^(r) = 5 l0 5{v) + dld^r)-\ (38) 
4(r) = -dld^r)' 1 . (39) 

We insert Eq. (^) into Eq. ( |56] ) and obtain the following integral representation 
of the vector potential: 

At (r , w) = no J d 3 s G\ m (r , s, u) j_ m (s,u), (40) 

where 

Gi m (r,s,u) = J d 3 s'6^r-s')G jm (s',s,u;) (41) 

is the (from the left) one-sided transverse Green tensor. Next we substitute in 
Eq. fl3"TD for the longitudinal 5 function the expression ( |3"9"D and find that 

0(r )W ) = ^ fd 3 s^± (42) 
— j 47r|r — s| 

The integral representation of the scalar potential can then be found, substituting 
in Eq. ( f4*2|) for the electric field the integral representation (^3). 

2.5.2 Commutation relations 

It is well known that A(r) and eoA(r) are canonically conjugated field variables. 
In order to calculate the commutation relation between them and the scalar 
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potential <f(r), we first note that A(r), A(r), and ^(r) are given by integrals of 
the type fll7|), but with A(r, uS), —iujA(r,uj), and <p(r,u), respectively, in place 
of E(r, uj). Using Eq. ( f40|) and following the lines outlined for calculating the 
commutation relations (ESI) and (130), we obtain 



and 



where 



|"/°°^5 G « (r ' r>) 

Tito J— oo C 



(43) 



(44) 



| d 3 s | dV 6t(r - s)G mn (s, a', u)6^{a' - r') 



(45) 



is the two-sided transverse Green tensor. Similarly, the commutation relation 
between the vector potential and the scalar potential can be given by 



£(r),i,(r') 



n 

vre 



8L / d 3 s 



du 



UJ 



G m )(s,r',a/) 



c 2 47r|r — s| 



where 



is the (from the right) one-sided transverse Green tensor, and finally 







0(r),Mr') 



(46) 



(47) 



(4* 



In order to further calculate the commutators fl43|) and fl46|), the Green function 
multiplied by the frequency must be integrated over frequency, which is quite 
similar to the commutation relation (|29T). 



3 Homogeneous dielectrics 



Let us first consider the electromagnetic field in an absorbing bulk material such 
that the permittivity can be assumed to be independent of space: e(r,u) = e(co>) 
= 6r(u;) + i€i(uj) for all r. In this case, the solution of Eq. (pi) which satisfies the 



boundary conditions at infinity is ]54] 



r ,u) 



q (u)g{\r 



(49) 
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where the notation q 2 (u) = (u 2 /c 2 )e(uj) is used, and 

iq(u)\r-r'\ 

g (\r-v'\,u;) = — -. (50) 

47r|r — r'| 

Substituting in Eq. ( ^3|) for the Green function the function given in Eq. ([49|) 
and integrating over frequency yields, together with Eq. (|I~8D, the Green function 
expansion of the operator of the electric field strength of the electromagnetic field 
in a dispersive and absorbing bulk dielectric. The Green function expansions of 
the magnetic field and the potentials can be obtained accordingly. 



3.1 Commutation relations 



We insert the Green function (|49D into Eq. fl29|) and find that 



h 

ne 



du — g{\r - r 

-oo C 



(51) 



[note that ekmjd^dj' (...) = 0, because of the antisymmetry of the Levi-Civita 
tensor]. We now substitute in Eq. (|5T"D for g(\r — r'|, uj) the Fourier expansion 



^(|r - y'\,u) = J 



d 3 k e ik ( r - r ') 
(2vr) 3 k 2 - q 2 (cu) 



(52) 



which enables us to calculate the u integral by means of contour integral tech- 
niques. Recalling the properties of e(Q) as a function of the complex frequency 
fl, we obtain (cf. 



, 41 



00 , u 
doj 



oo c 2 k 2 — q 2 (uj) 



(53) 



Since the result is independent of k, the k integral simply gives a 5 function, and 
we eventually arrive at 



£i(r),4(r') 



ih 

eo 



e ikm d r m 5{r - r'). 



(54) 



From Eqs. ([54]) and ( p0|) we see that the quantization scheme yields exactly the 
equal-time electromagnetic-field commutation relations that are well established 
in QED. 

In order to calculate the commutator [A(r), Aj(r')], we combine Eqs. (f4"3|), 
(flop, and (|49|). After a partial integration the derivatives in Eq. (^9|) give rise to 
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a divergence of the transverse S function in Eq. ([43|) which vanishes, d'-S^r) = 0, 
so that 



— / d 3 s I <Ps'd- (r-s) 

7Te 



/°° CO 
duj— g(\s- s'\,u) 5 mj (s' - r'). 
-oo C 



(55) 



We again make use of Eqs. ( |52"D and (|55D, and after some straightforward calcu- 
lation we obtain 



A(r),^(rO 

Quite similarly, it can be proved that 

0(r),Aj(r' 



0. 



(56) 



(57) 



Again, the correct QED equal-time commutation relations are preserved. Note 
that from Eqs. fl32|), (pi), and (|7|) it follows that 



-^(r-r). 



(5? 



3.2 Relation to earlier work 

To make contact with earlier work, we first note that, according to the Helmholtz 
theorem, the noise current j can be decomposed in a unique way into a transverse 
and a longitudinal part, 



= j (r,w) +j (r,w) 



where 



*-L(ll), x 



d 3 s^ (ll) (r-s)f(s,o;). 



— j 



(59) 



(60) 



Since the permittivity is assumed to be independent of space, the decomposition 
of j corresponds to a decomposition of the bosonic (basic) field f into a transverse 

and a longitudinal part. We use Eq. (|18D and can relate j ^ to f~ L< -"' ) as 



?-L(ll), v 

J [r,w) 



fi c z V Treo 



where 



(61) 



(62) 
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The commutation relations ([19]) and (|20| ) obviously imply that 

= 4 (ll) (r-r')%;-u/), 



(63) 



and 



/^(r,u;),(jf(rV)^ 



0. 



(64) 



(65) 



Obviously, the motion of the transverse and longitudinal degrees of freedom are 
governed by their own Hamiltonians, as can be readily seen from the Hamiltonian 



(21), which can be rewritten as 



H 



(66) 



with 



d r / duTiuj 



fJ-(ll), 



.fJ-(ll)( 



(67) 



From Eqs. (f40|), (|4l|), and (f49|) it can be seen that, after partial integration, 
the derivatives in Eq. (^) do not contribute to the vector potential, because of 
the vanishing divergence of the transverse 5 function. Therefore we may write 



(68) 



Ai(T,u) = Ho J d 3 s 

x J dh' 5i j (r-s)g(\s-s'lu)l j (s\uj) 

or, changing the variable s as r' = r — s + s', 

A*(t,w)=Ho / d 3 r' g(\r-r'\,u)j. (t',uj), 



(69) 



which is nothing but the representation of the transverse vector potential given 

in |^D[ (if we identify /ioj with j in |4"0|]). Similarly, from Eqs. (^) and ([4"2]) 
and the Green function (]49D [together with Eq. (|50D1 it can be derived that 



47reoe(u;) 



3 P( s -^) 



r — s 



d 6 s 



(70) 
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where 



p(r,uj) = (iu) V ■ j (r,w) 

= -^J- f ,HV-f"(r, W ) (71) 

[cf. Eq. (|TB|)]. Note that in the case under study from the commutation relations 
fl65D and Eqs. (|69D and ( [f0| ) [together with Eq. ([71])] it is immediately seen that 
vector and scalar potentials are commuting quantities, i.e., the commutation 
relation ( |57D is fulfilled. 



4 Dielectric Interface 



4.1 The Green function 

The determination of the dyadic Green function for three-dimensional configura- 
tions of dielectric bodies is a very involved problem in general, and only for rather 
simple configurations the Green function has been calculated explicitly. Such a 
configuration, which can be thought of as being the basic element of multilayer 
dielectric structures, is composed of two infinite half-spaces (Vi and V 2 ) with a 
common planar interface such that 



6x{UJ) 
62 M 



if 
if 



r G Vi, i.e., z < 0, 
r G V 2 , i.e., z > 0. 



(72) 



Following |^0|, ^8], |54| , we write the solution of Eq. ( p4|) in the form 

(G%(ry,u)+B%(ry,u) if r,r'GV a , 
\T t f(r, r >,u) if reV a , v'eV al (a^a') 



(73) 



a' = l, 2), where Gfj(r, r', uj) is the Green function ( f4~9|) [together with Eq. (|50|)1 



for the bulk material with e a (uj), 



G%{vy,u) 



dm + 8 ijq l{u) q~ 2 {uj)g a {\v 



with 



,0J 



e ig a (a»)|r-r'| 

4/rlr — r'l 



(74) 



(75) 



[q^{ui) = (lu 2 /c 2 )e a (u)]- The functions R?Ar,r',u>) and T°j a '(r, r', uj) describe the 
effects of reflection and transmission at the interface and obey the homogeneous 
equations 



d r d r - 5- 

u i u m U V 



v. 







(76) 
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(r, r' G V a ) and 



d r d T -5- (A r + a 2 



rpctOL 

mj 



r,r ,u) 







(77) 



(r G V a , r' G V a >, a / a') together with the boundary conditions such that they 
vanish at infinity and the tangential components of the E and the H field are 
continuous at the surface of discontinuity In order to determine the functions 
i?^(r,r' ,u) and T°- a (r, r', u), it is convenient to take advantage of the transla- 
tional invariance of the system in the (x, y) plane, which enables us to expand 
the Green function (and G"-, i?"-, and T°j a ) as follows: 



d 2 kn 



ik n-( r i[- r ;) 



(7i 



where kn = (k x , k y , 0) and ry = (x, y, 0) are two-dimensional vectors in the (x, y) 
plane. Using the general formulas for multilayer structures given in [j54] (see 
also p)Tf), after some manipulations we obtain, for the reflection functions R\ 



19. 

i] 



RC* = J_ e iM\z\ + \z'\) 



2(3 a 




_l_ r s _y_ 

~ act' t.2 



(79) 



pa 

^xy 



2(3 a 



s i^„(|*|+|*'|) 



x 



k k 

2 Kg; fry 

~ Pa kl 



k k 

s to x rvy 



RZ 



A e ^(W + k'l)^i[_ sign( ^ )/3afe]) 



pa pa 

^yx xyi 



Ryy — Rxx(kx ky)i Ryz ~ ^xzi^x *~* k y ) , 



pa 



pa 

1 h v ~ 



-R1 



yz 

pa 
^zy 



pa 



1 if3 a (\z\ + \z'\) r aa' h 2 

2/?a ql 11 ' 



and accordingly for the transmission functions T^ a = Tf^ a ' (k\\, cu; z, z') 



52) 
i3) 

15) 



iaa' __ % i/3 a \z\+i/3„,\z'\ 
xx " 2(3 a 



it 



h . ' 



kl 



x 



q a q 



(86) 
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rnaa! _ _J_ i(3 a \z\+i/3 a ,\z'\ 

xy 2(3 a 

t p k h h h 

b aa' g a ^x^y , s l^x^y i / Q7 \ 

— — PaPa'—pr l aa'—j^T I ' \° ') 

q a q a > in tii 

T«? = -!-e iMzl+i ^' lz ' l ^^sign(z')p a k x , (88) 
2 p a q a qa' 

T™' = T™', (89) 
Ttf = T™'(k x ~ k y ), T y a ? = T:?(k x ~k y ), (90) 

T^' = ^e iMzl+l ^' lz ' l ^sign(z')p a ,k x , (91) 
2 p a q a q a ' 

Tg = T™'(k x ~k y ), (92) 

J^' _ * ^ffaM+ifl,,!*'! V ft2 (-g 3 ^ 
2/3 Q 9a9a' 

where a' = 1(2) for a = 2(1), and 

q a = q a (u), 

Pa = A» = y/ql(u) - kp Re Pa > 0, Im /3 Q > 0, (94) 

with r aa , = r' a ,(u;) and £* a , = t aa ,(u), q = p,s, being the generalized reflection 
and transmission coefficients for the p- and s-polarized components of the elec- 
tromagnetic field, which are defined by 

q _ Pa - laa'Pa' _ q /qr\ 

7^ = —, 7^ = 1, (96) 
tL' = \flL>(l+r 9 aa ,) = § L t q a , a , (97) 

Pa' 



(for details, see ]5~|). 



It should be pointed out that - in contrast to the usually considered external 
current - the noise current j as given in Eq. flTE| ) jumps at the interface, i.e., at 
z = 0, which obviously implies the existence of a surface noise charge density, 
and hence the normal component of eE is not continuous at z = 0. For more 
details and a derivation of the fields E and B by direct solution of the Maxwell 
equations, the reader is referred to App. |B[ 
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4.2 Commutation Relations 



As shown in App. 0, the functions R?j(r, r', u) and T^ Q '(r, r', u) [Eqs. d? 1 
obey the relations 



where 



duj—R?Ar,r',u) = d r s d r , R a (r, r') 



—oo C 



.2 l J 1 



oo C 



.2 V 



R a (r,r') 



d 2 ky 
(2^ 



ik|i-(i-||— r') 



x / ^ J_ e ^<M+i*i)I«s: 

3 C J 2/? Q <? 2 



(98) 
(99) 



(100) 



rpaa> 



r, r 



d 2 k„ 



-, ik ir( r ii- r !i) 



(2vr) 2 

, a; i 
aa; 



c 2 2/3 Q 



3 i/3 a |z|+i/3 Q ,|2;'| 



(101) 



Further it can be shown that 

r^Gg(r,r>) 

J —oo C 

= i7r<%5(r - r') - dld r -G a {r, r'), 
where we have used the Weyl expansion of g a (\r — r'\,uj) |54| to obtain 
G a {r,r') 



(102) 



c 2 2/5 a g 2 



(103) 



When we substitute in Eq. (|29| ) for (jy(r, r', a;) the actual Green function (|73|), 
take advantage of the relations (0), fl99"D , and ( |102| ), and recall that €k m jd r m dj (...) 
= 0, then we can readily prove that the fundamental QED commutation relation 



^(r),4(r')l = ~-e ikm d r m 5(r - r') 



(104) 



is preserved. 

In order to prove that the commutations relations between the potentials are 
correct, we first write Eqs. (|43|) and fl46"D [together with Eqs. (^) and (|47D 1 as 



Mr), Mr') 
y3(r),i,(r') 



(105) 
(106) 
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where 



4 ) (r,r') = ^/rf 3 s/rfV X S(r-s) 

<ko—G mn {s, s', u)5 nj (sf - r'), 

-oo C 



with 



and 



(107) 
(108) 

, , , • ( 109 ) 
47r|r-s| v ' 

We use the Green function ( |73|) and take again advantage of Eqs. Q98|) , (p9|), and 
( |i02D to derive 

2 



Xim(r-s) =5^(r-s) 

(2)/ \ &im 



r (fe)/ /\ ^ I 

!„■„■ Mr, r J = 

^0 a 



E £ ^xS(r-s) 



X 



iTiS mn 5(s - s') 



(-G a (s,s , )+fi a (s,s'; 



ih 



- E / - " O + r ') 

e a=l JVa 



ih 

e 



y ^ 3 sxSS( r -s)^-(s-r') 



(110) 



[a' = 2(1) for a = 1(2)]. The last line follows from the line before last, because of 



i$\ry) = — E/ <* 3 sx&v- £ 

7re V J v a 



n 



/ dV fc<9*' (-G a (s, s') +iT(s, s' 



+ / dV 



d s m d s n T™ (s,s') 



^(s'-r') 



(111) 



(App. 0). Recalling the definitions (|108Q and ( |109| ) of xill an d Xim? respectively, 
from Eq. ( |110[ ) it is easily seen that 

= ^ J d 3 s5j- m (r - s)5^(s - r') 



iff 



fo 



5h(r-r') 



(112) 



17 



and 



d r I {2) (r r') = - [ d 3 sd r ^ (S ~ r/) 
a x i l3 [v,v) e J aso m 47r|r _ g| 



o, 



(113) 



i.e., 



and 



Air), Mr') =^(r-r') 



e 



(114) 



^(r),i,(r) 



0. 



(115) 



Hence the theory yields the correct equal-time commutation relations for both 
the fields and the potentials. 



5 Conclusions 

We have developed a quantization scheme for the electromagnetic field in a spa- 
tially varying three-dimensional linear dielectric which gives rise to both disper- 
sion and absorption. Based on the classical phenomenological Maxwell equa- 
tions, the dielectric is described in terms of a complex frequency- and space- 
dependent permittivity, which satisfies the Kramers- Kronig relations, and fluctu- 
ating current and charge densities are introduced in order to be consistent with 
the dissipation-fluctuation theorem. The noise current and charge densities can 
be thought of as arising from an additional noise polarization in the constitutive 
equation between the dielectric displacement vector and the vector of the electric 
field strength. The resulting inhomogeneous Maxwell equations are then trans- 
ferred to quantum theory, and the noise polarization is specified such that the 
fundamental equal-time commutations relations of QED are preserved. 

From the inhomogeneous Maxwell equations together with the boundary con- 
ditions at infinity it follows that the electromagnetic field operators can be related, 
through the dyadic Green function of the classical problem, to a bosonic field that 
represents the elementary (energy) excitations of the overall system. This integral 
representation can be regarded as a natural extension of the familiar mode expan- 
sion in free space or in cavity-like systems with perfectly reflecting walls. Vector 
and scalar potentials are introduced in the usual way, and their integral repre- 
sentations are derived, which can be used in order to couple the electromagnetic 
field to additional atomic sources embedded in the medium. The fundamental 
equal-time commutation relations are studied, and it is found that some of them 
can be calculated without kowledge of the explicit form of the Green function 
in order to prove the consistence of the quantization scheme with QED. Others 
require this knowledge, because single frequency integrals of the Green function 
remain to be calculated. 
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The determination of the Green function is - similar to the determination 
of the mode functions in a mode-expansion approach - a purely classical prob- 
lem. Its solution is very difficult in general, and only for simple dielectric-body 
configurations the Green function has been calculated so far. For a homoge- 
neous dielectric the Green function is well known. We have used it and explicitly 
shown that the quantization scheme outlined here yields exactly the fundamental 
equal-time commutation relations of QED. We have further shown that earlier 
results derived in [[!(], [1(J for the transverse part of the electromagnetic field are 
contained in our theory. 

An example of an inhomogeneous medium for which the Green function is 
known is a configuration of two infinitely extended dielectric bodies with a com- 
mon planar interface. We have also used this Green function and explicitly proved 
that the quantization scheme is consistent with QED. Again, recent results given 
in [[H], fll]] for paraxial light propagation are recognized. It is worth noting that 
the contributions to the Green function which result from the reflections and 
transmissions at the interface do not contribute to the equal-time commutation 
relations. Since this is expected to be true also for more complicated configu- 
rations, such as multilayer structures, the equal-time commutation relations are 
expected to be preserved also for these configurations. 

The quantization scheme developed in this paper can be regarded as the 
basis for studying the interaction between radiation and atomic systems in the 
presence of three-dimensional configurations of dielectric bodies with dispersion 
and absorption. In this case, additional atomic sources must be introduced into 
the theory, and the coupled equations of motion for the atomic variables and the 
electromagnetic-field variables in the Green function expansion (similar to the 
photonic variables in a mode expansion) must be tried to be solved. 
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We multiply Eq. (|116|) from the right by G*As,r' ,u) and integrate over s. We 
derive, on integrating by parts and taking into account that the Green function 



A Proof of Eq. (|27|) 




(116) 
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vanishes at infinity, 



— J d 3 se(s,u})G H (s,r,u)G u { j (s,r',Lj) 

= J d 3 s[(dfd s n -5 ln A s )G ni (s,v,uj)]G^(s,v',aj) 



-G*ij{r,r',u}). 



-J d 3 s [d s n G ni (s,r,u)] [dfG^sy^ 

+ Jd 3 s [dlGuis^u)] [^(s/.w)' 
-G* tj (r,r',uj). 



(117) 



Now we take the complex conjugate of Eq. ( |117| ), make the interchanges r <-> r' 
and % j, and subtract the resulting equation from Eq. ( 117 ). In this way we 
arrive at 



J d 3 s ei(s, lo)G u (s, r, w)Gy(s, r', u) 



2i 



Gu{v',v,u)-GUv,v\lo) 



Recalling the symmetry property |p2 | 

Gji(r',r,uj) = Gij(r,r',uj) 



(118) 



(119) 



and combining Eqs. (|118j) and (119) then yields Eq. 



B Direct solution of the Maxwell equations 



In order to directly solve the Maxwell equations (|10|) - (|T^) together with the 
permittivity in Eq. ([72]), i.e., without using the Green function, we expand E(r, to) 

as 

= / ^LE(k |h o;;^ e *ir r il (120) 

[cf. Eq. (|78j ) and the comment made there]. Obviously, similar expressions hold 
for B and j. Choosing a reference system in which ky is parallel to the x axis, 
the Maxwell equations (|12"D and (|13|) then yield 

-d z E y = *ujR x , (121) 
d z E x -ik\\E z = iuB y , (122) 
ik\\E y = ituB z , (123) 
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IU1 

-d z B y = --eE x + mx , (124) 
= -^E y + fi 3_ y , (125) 

= -^ek z + ^j_ z , (126) 

where, for notational convenience, we have omitted the arguments kn, uj, and z 
of the fields E(k\\,uj; z), B(k[|, uj; z), and j(k||, uj; z), and we have used the same 
notation for the rotated reference system as for the original one. According to the 
Maxwell equations (|l"2"D and flT3|), the tangential components of E and H must 



be continuous at the surface of discontinuity ]58| . 



From Eqs. ( |121| ), ( |123| ), and ( |125| ) together with the boundary conditions we 
find that E_ y obeys the equation 

d 2 z E y + (^e - k^j E y = -iu^[ y7 (127) 

and E_ y and d z E_ y are continuous at the plane z = 0. B^. and B_ z can then be 
obtained from E y , using Eqs. O and (|T23|) . Similarly, from Eqs. ( |i~22l) , (|T2l) , 
and ( |126| ) together with the boundary conditions it follows that B_ y satisfies 

d z R y + l^je - kj] B y = ik^l - fi dj x , (128) 



and B_ y and e 1 (d z B_ y + fi j x ) are continuous at the plane z = 0. Knowing B_ y , we 

can find and E z from Eqs. Q) and (P5|) . 

Thus, Eqs. (|121|) - ( |126| ) together with the boundary conditions at infinity 



and the conditions of continuity of the tangential components of E and H at the 



interface provide us with unique solutions for E and B. Solving Eqs. ( |127| ) and 
128|) under the conditions mentioned (and going back to the original reference 



system) yields the fields E and B in full agreement with Eqs. (|23|) and (^) and 
the Green function from Sec. |4.1| . It should be emphasized that, as can be seen 
from Eq. ( |126| ), the normal component of eE is not continuous at the plane z — 
0, because B_ y is continuous and j makes a jump at z = [see Eq. (|T8| ) together 
with Eq. (|72D]. 



In the derivation of the Green function in [51] it is stated that the continuity 
of the normal component of eE at the plane z = is chosen as a boundary 
condition. However, one can verify that it is the continuity of eG Z i, i = x,y, z, that 
actually matters, and this can be deduced from the continuity of the tangential 
components of E and H. From the continuity of eG z i it does not follow that eE_ z 
is continuous, though the reverse is true. For example, when (as in our case) 
eG zz contains a term proportional to 5(z — z'), which is continuous at the plane 
z = 0, then eE_ z jumps there, because it contains a term proportional to j , which 
makes a jump at z = [see Eq. (pif) together with Eq. ([72])]. 
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C Proof Of Eqs. (|9§) and (|99D 



In order to prove Eqs. (|9~8D and (|9~9]) for Rfj(r, r', to) and T^ a '(r, r', a;), it is suffi- 
cient to prove them for the Fourier components i?^- (ky , a;; z, 2;') and T^ a (ky , cu; z, z') . 
From Eqs. (|79|) - (p5|) and Eqs. (|86|) - (|93|) , respectively, the functions i?°-(ky, w; 2, z') 
and T^ a '(ky,u;; z, z') are seen to consist, in general, of two parts. One part is as- 
sociated with the p-polarized and the other part is associated with the s-polarized 
electromagnetic field. Using contour integral techniques (cf. [0, 0), it can be 
proved that the latter does not contribute to the integrals on the left-hand sides in 
Eqs. (|9~8|) and (|9~9"D , so that we are left with the contributions from the p-polarized 
field only. 

Let us consider, e.g., R^ x (k^,u; z, z'), Eq. (^). From 



2f3 a e ql V Nk 



-00 C Zp a q a K|| 

^ du ^ 'iM\AMA) r lf k l (129) 



c 2 2(3 a ql 



Eq. ( P5| ) can readily be proved correct for R® x (r, r', u). In exactly the same way, 
one can show that Eq. (j98D also holds for the remaining functions -R^ (r, r', cj). In 
order to show that Eq. (^) is valid, it is helpful to perform the integration over 
u in Eq. ( |101|) explicitly. Using the relation 



^■/"s^jR' (130) 



we have 



00 



c 2 2/3, 




'no' 



00 c 2 2(3 Q (3l-k 2 (3l,-k' 2 \ q a q 



(131) 



Since the equation j3 2 — k 2 = has no solutions in the upper complex frequency 
half-plane [35, 41], the integrand in the u integral in Eq. ( |131|) has no poles there. 



Treating the single pole at to = as principal value, we obtain 

3 c zp a \ q a q a ' 
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00 dk e ikz r°° dk' e ik ' z ' 



, 2tt A; 2 + A; 2 7^oo 2tt A; 2 + k' 2 

x 2 OT A;|| | ; ^ ,(0) 1 . (132) 

Let us now turn to the left-hand side in Eq. (|99D and consider, e.g., T""' (ky , cu; z, z') 
Eq. (H). We derive 



oo 



duJ W 1 e i(3 a \z\+i[3 a ,\z>\ Kaj PaAs 
oo C 2 2f3 a <la(la> k 2 

~dk eitz r»«>, ~ „ 

-co 27T J —oo 27T J —oo C 

i 2i(3 a 2i(5 a < t p aa , [3 a f3 a > 
X 2(5 a (3l-k^l,~k' 2 q a q al k 2 



3 2tt A| + A; 2 7-oo 2vr A; 2 + A;' 2 
x2zttA;|| f . I , (133) 

which exactly agrees with Eq. ( |132| ), i.e., Eq. fl99|) is proved correct for T"^'(r, r', u;). 
Eq. (|132|) for the other functions T°j a (r, r', u) can be given in the same way. 



D Proof of Eq. ( |TTTl ) 



We perform in Eq. ( |1 1 1| ) a partial integration over s' and obtain, on recalling 
that ^'<5i(s , -r / ) = 0, 



'n "raj 1 



4?(r,r') 



7re 



Z! / ^ 3 sxS( r -s) ds' x ds' y 



xd s m [-G a (s, s') + R a (s,s') 
-T^'(s,s')]5i(s'-r') 

Using Eqs. (|100|), ( TO , and ( |T03|) , we derive 

-G°(s,s') + /T(s,s') -f Qa '(s,s') 



s' =0 



/ 



d 2 k|| 

(2^ 



ik||-(s||-B{,) 



s i=0 

00 ^ u i 

— oo C 2 2/^q, 



(134) 
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1 rL, e 



rf 2 kn ik,,.(8||-B',) r°° . w z 



(2tt) 2 7-oo c 2 2/3, 



V 9a 9a Pa' V 6 « 9«<?a 

2 




« kn ikN-fsM-sn f°° , uj i 



(2tt) 2 7-oo c 2 2/3 a 

"^hl)^ ' (135) 



To obtain the second equation from the first one, we have used Eq. (|97|) for 
t p aa ,. The third equation gives zero because of the integration over u. Combining 
Eqs. O and Eq. O yields Eq. ([TTTD . 



*Permanent adress: Institute of Physics, National Center for Sciences and Tech- 
nology, 1 Mac Dinh Chi St., Dist. 1, Ho Chi Minh city, Vietnam. 
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